Nonequilibrium Relaxations and Aging Effects in a Two-Dimensional Coulomb Glass 
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The relaxations of conductivity have been studied in the glassy regime of a strongly disordered 
two-dimensional electron system in Si after a temporary change of carrier density during the waiting 
time tw ■ Two types of response have been observed: a) monotonic, where relaxations exhibit aging, 
I. e. dependence on history, determined by t^ and temperature; b) nonmonotonic, where a memory of 
the sample history is lost. The conditions that separate the two regimes have been also determined. 

PACS numbers: 71.55.Jv, 71.27.+a, 71.30.+h 



The dynamics of glasses and other systems out of equi- 
librium is one of the most challenging and rapidly evolv- 
ing topics in condensed matter research Although 
glassy behavior may dominate the low-temperature prop- 
erties of many complex materials near quantum phase 
transitions such as the metal-insulator transition 
(MIT), quantum glasses are less understood than their 
classical counterparts. Experimental studies of charge 
or Coulomb glasses Q , which are of particular relevance 
to the MIT, have been scarce 0, 0, @|, and limited to 
insulating systems far from the MIT. Recent observa- 
tions 0, di @| of glassiness in a two-dimensional electron 
system (2DES) in Si open up opportunities for exploring 
glassy phenomena over a wide range of all the relevant 
parameters. The carrier density n^, for example, can be 
varied continuously on the same sample from the deep 
metallic into the deep insulatingregime. The onset of 
glassiness was, in fact, found 0,0, 0] to take place at a 
density Ug > ric {ric ~ the critical density for the MIT), 
i.e. on the metallic side of the MIT, consistent with the- 
ory 0. 

The manifestations of glassiness in a 2DES for rig < rin, 
as demonstrated by resistance noise measurements 0, Q j 
include a dramatic slowing down of the electron dynamics 
and correlated statistics consistent with the hierarchical 
picture of glasses. Furthermore, studies of the relaxations 
of the conductivity a P , following a rapid change of Ug , 



show that the equilibration time Teq 11| diverges expo- 
nentially as temperature T ^ 0, suggesting a glass tran- 
However, even at T that are not too low 
eq exceeds easily not only the experimen- 



sition at Tg = 



{e.g. - 1 K), r, 

tal time window but also the age of the Universe. This 
makes it relatively easy to study the out-of-equilibrium 

relaxation of a at times t ^r^q, where one expects to find nc(10^^cm~^) fa 4.5 [9| . The samples and the standard ac 



The most puzzling feature of the relaxation in Ref. [9|, 
however, is that the system equilibrates only after it first 
goes farther away from equilibrium. This "overshooting" 
(OS) of (To is discussed in more detail below. 

A key characteristic of relaxing glassy systems is the 
loss of time translation invariance, reflected in aging ef- 
fects ^M,- The system is said to exhibit aging if its 
response to an external excitation depends on the system 
history in addition to the time t. Some history depen- 
dence was observed in a 2DES , but there have been no 
systematic and detailed studies. Here we report the first 
such study, where, for a given rig, (j(t) was measured after 
the temporary change of Us during the waiting time tw 
The measurement history was varied systematically by 
changing and T. Two types of response have been ob- 
served: a) monotonic, for relatively "small" excitations, 
where relaxations depend on t^, i.e. the 2DES exhibits 
aging; b) nonmonotonic, for sufficiently "large" excita- 
tions, where a{t) overshoots ctq and relaxations no longer 
depend on t^ (memory loss). The conditions that sepa- 
rate the two regimes have been identified and discussed. 

Measurements were carried out on a 2DES in Si metal- 
oxide-semiconductor field-effect transistors with a 50 nm 
oxide thickness, and a relatively large amount of disor- 
der (the 4.2 K peak mobility « 0.06 m^/Vs with the 
substrate (back-gate) bias [3] Vsub = —2 V). The sample 
length X width were 1x90 htd? (sample A) and 2 x 50 /um^ 
(sample B). The same devices were used in Ref. Q, and 
they were also very similar to the samples used in noise 
studies . The density is controlled by the gate volt- 
age Vg, such that ns(10"cm-2) = 4.31(Fg[V] - Vt), 
where the threshold voltage for forming the 2D layer 
Vt = 6.3 V for Knf,_= -2 V; ng(10"cm-2) w 7.5 and 



properties common to other types of glasses. Indeed, for 
t < T, 



pes o 

eg, the relaxations obey [9] (T{t,T,ns)/a'o{T,ns) oc 



^-a(n,) exp[-(t/r(r, n^))'3("=)] with r diverging as T ^ 
{ao{T,ns) - equilibrium a at a, given T and rig; < 
a{ns) < 0.4, 0.2 < l3{ns) < 0.45). Such scaling has been 
observed in spin glasses above Tg The dependence 

T oc exp(7ny^) (7 - a proportionality constant) is a 
strong indication for the relevance of Coulomb interac- 
tions between 2D electrons [is'] in the < Ug regime. 



lock-in technique (typically 13 Hz) that was used to mea- 
sure a have been described in more detail elsewhere 0. 

The experimental protocol [Figs.[lja) and[ljb)] starts 
with the 2DES, induced by an applied gate voltage Vq, 
in equilibrium at 10 K. The sample is then cooled to the 
measurement temperature T. Unlike glasses with Tg ^ 0, 
this does not result in any visible relaxations on our ex- 
perimental time scales, and ao{Vo,T) is thus estabhshed. 
Vg is then switched rapidly (within 1 s) to a different 
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FIG. 1: (Color online) (a) Sample A. a{t) for Vo = 7.2 V 
K(10"cm-2) = 3.88 < nc], Vi = 11 V K(10"cm-2) = 
20.26], t» = 1000 s, and T = 3.0 K. (b) The corresponding 
experimental protocol: Vg{t) and T{t). The results do not 
depend on the cooling time (varied from 30 minutes to 10 
hours), nor on ti and t2 (varied from 5 minutes to 8 hours 
each), (c) Sample B. a{t) for Vo = 8 V [n^lO^^cm'^) = 
7.33 ^ rig], Vi = 11 V, — 1000 s, and several T, as shown. 
T = 2.3, 1.8, 1.3, 1.0, 0.75, 0.5 K data were omitted for clarity. 



value Vi, where it is kept for a time ty^. Finally, it is 
changed back (within 1 s) to Vb, and the slowly evolving 
a{t,Vo,T) is measured. We note that t = is defined 
as the time when Vg reattains its original value Vq. An 
analogous protocol has been used often in studies of other 
glassy materials [H|. In general, we expect the larger val- 
ues of or T to correspond to stronger excitations of 
the 2DES, because the system will have more time or 
more thermal energy, respectively, to wander away from 
its original equilibrium state during the time 

Figure [Ijc) shows some relaxations measured for a 
fixed at different T for a given Vq and Vi. The data 
are not shown for the first few seconds (comparable to 
our sampling time ^ 1 s), when a changes rapidly, con- 
sistent with the (Jo{Vg) dependence In general, the 
observed behavior is very complicated, but it can be di- 
vided roughly into two types. At higher T, a{t) first 
goes below cto , and then even farther away from equilib- 
rium. Eventually, it changes "direction" and approaches 
ctq. Similar nonmonotonic behavior was observed after 
exciting 2DES out of equilibrium by a rapid change of 
Us 0. At lower T, on the other hand, there is no OS of 
equilibrium and a approaches ctq monotonically, similar 
to observations in other electron glasses 5]. 

a(t) were also measured for different at a fixed T 
and given Vq and Vi . Some typical data, including those 
for two different Vo, are shown in Fig. [51 In general, the 
waiting time tw has a dramatic effect on the relaxations, 
both qualitatively (OS vs. no OS) and quantitatively 
{e.g. more than two orders of magnitude difference in 
a{t) in Fig. EJd)), including instances of OS when t^- 
dependence vanishes [Fig. [^Jb)]. a{t) were also investi- 
gated after all the electrons have been removed from the 



2D layer during the time t^,. In a 2DES in Si, unlike other 
electron glasses, this may be accomplished simply by re- 
ducing Vg by a sufficient amount. If the slow dynamics 
is dominated by intrinsically glassy behavior of the inter- 
acting electrons rather than by their response to extrinsic 
slow degrees of freedom {e.g. motion of nearb y im puri- 
ties, although exponentially unlikely at our T til), one 
expects that <7{t), measured after the electrons are rein- 
troduced with the initial Vg — Vq, will not depend on tw 
The measurements confirm this expectation [Fig. [2l^f)], 
providing further support for the glassiness in the 2DES. 

A detailed analysis of cr(t) in the case of OS [e.g. high- 
est T curves in Fig. [^c)] finds the behavior very similar 
to that reported in Ref. 0]. Therefore, we do not show 
the fits to the data, but instead summarize the main 
results. For example, at t just before the minimum in 
a{t) a{t,T^s)/ao{T,n,) oc exp[-(VTjr, n,))'^^"^)], 
with T and /3 [18| similar to those in Ref. [9|. After the 
minimum in cr(t), the approach to equilibrium at the 
longest t is described by a simple exponential process 
a/ao = [l-A{T) exp(-t/Tp)] with 1/Teg oc exp(-S^/r) 
{Ea « 57 K as before [9|, [ill). The main difference be- 
tween the equilibrations that follow the two experimen- 
tal protocols is that the prefactor A did not depend on 
T in Ref. whereas here A{T) decreases with decreas- 
ing T and vanishes at a finite T (Fig. [3Ka) inset). This, 
at least partly, reflects the observation that, in Ref. 
the minima in the a{t) curves only move to lower val- 
ues of a and to later t with decreasing T, whereas in 
the current protocol, the different T curves at t after the 
minimum [Fig. [ijc)] are no longer parallel to each other, 
the minima become "shallower" , and eventually vanish 
at a finite T. The latter is shown in Fig. [3Ja), where 
the T-dependences of the amplitude of the relaxations 
at the minima, defined as [log{a / ao)]mim are plotted for 
different i^,. The amplitudes extrapolate to zero at T 
that are in a good agreement with those obtained from 
A{T). Thus the data seem to indicate that here the OS 
does not simply fall out of our experimental time win- 
dow at low T as in Ref. [§| but rather that it vanishes 
at a finite T. The a{t) curves at lower T support this 
conclusion: their monotonic approach to ctq [from above 
in Fig. [ijc)] is consistent with a power-law form at the 
shortest t (lowest T) and, at the longest t, with a simple 
exponential process with the same Ea « 57 K. 

The procedure shown in Fig. [Slja) was applied to all 
the data in order to determine the values of (T, tyj) where 
OS disappears. The results [Fig. Eljb)] clearly divide the 
phase space in two, such that the OS is not observed 
for smaller values of (T, i^) {e.g. at 1 K, not even at 
t ~ 2 X 10^ s). This confirms our earlier speculation 
that the OS occurs for sufficiently strong perturbations 
out of equilibrium. We note that the relevant time (en- 
ergy) scale is given by Teq{T) [Fig. EKb)]. In order to 
understand this result better, we discuss the experiment 
in more detail. During t^, after Vg is switched rapidly 
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FIG. 2: (Color online) a{t) for different t„ at fixed T, as shown, (a)-(b) Sample B; Vo = 8 V, Vi = 11 V; (c)-(e) Sample B; 
Vo = 7.2 V, Vi = 11 V. (f) Sample A; Vo = 7.2 V and Vi = 6 Y{ns = 0) < Vr = 6.3 V, i.e. all electrons are removed from the 
2D layer during tw 
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FIG. 3: (Color online) (a) log(cr/(To) at the a{t) minima (data 
with overshooting) vs. T for different t^. Sample B; Vo ~ 
8 V, Vi = 11 V. The lines guide the eye. Inset: A{T) for 
tu, = 1000 s. The solid line is a linear fit. (b) The symbols 
(solid and open from devices B and A, respectively) show the 
values of (T, t^,) where the overshooting vanishes. Red dots 
and green triangles correspond to Vb = 8 V, Vi = 11 V, and 
Vo = 7.2 V, Vi — 11 V, respectively. The blue dot-dashed line 
guides the eye. The purple dashed line represents Teq{T). 



from Vo to Vi, a{t^ Vi,T) evolves with t as described in 
Ref. . Some small amplitude relaxations are still visible 
[albeit not on the scale of Fig. [TJa)] even at the highest 
Vi = 11 V where Ug > Ug, but ctq is still so small that 
kpl < 1 [kp ~ Fermi wave vector, I - mean free path) [l9| . 



If the 2DES equilibrates at a new state corresponding to 
Vi, i.e. if Teq <C (jeq docs not depend on 0]), then 
the rest of this protocol, namely the switching of Vi back 
to Vq and the subsequent relaxation, are equivalent to 
the protocol employed in Ref. Q. Hence both here and 
in Ref. [§] , the OS is observed only when the 2DES is ex- 
cited out of thermal equilibrium. Of course, if Teq <C tu,, 
then any further increase in will have no effect on a{t), 
as indeed observed [Fig. [2l^b)]. On the other hand, when 
tw ^ Teq [monotonic response regime in Fig. Eljb)], <7{t) 
clearly depend on history (aging effect), i.e. on t-u, [see 
e.g. Fig.[2)Ic)] during which the system relaxes away from 
its initial equilibrium state determined by Vq and towards 
a new equilibrium state determined by Vi . It can be also 
said that the system has a memory of the time it spent 
with Vg — V\. This is very similar to spin glasses, for 
example, where T or a magnetic field play a role of Ua in 
our experiment. Finally, in the crossover regime tw lS Teq 
[a region between the two dashed lines in Fig.[3fb)], both 
some OS and tu,-dependence are observed [Fig. [Hd)] . 

The nonmonotonic response of cr, which occurs when 
the system is taken out of thermal equilibrium, has been 
observed for all applied 1 <\AVg\^\Vi ~ Vo\ (V)< 4.7 
(I Aus I /tIs « 20 — 100%). Since, by reducing Ug, the 
2DES goes from a metallic regime, with no metastable 
states, to a glassy, and then glassy insulating regime, it 
seems reasonable to assume that the density of metastable 
states increases as Ug is reduced by a large amount, re- 
flecting the increasing complexity of the free energy land- 
scape. The 2DES are also known to be relatively poorly 
coupled thermally to the surrounding lattice, especially 
in Si, and thermal equilibrium at very low T (below 
'--^l K) is achieved mainly by conduction through the con- 
tacts and the leads. This and the available data suggest 
that the rapid change of Vg is an approximately adiabatic 
process. In that case, the reduction of Hs, accompanied 
by an increase in the density of metastable states, would 
lead to the cooling of the 2DES. As a result, a would de- 
crease with t {da/dT > in the regime of interest) below 
(To(T). This initial cooling of the 2DES and the eventual 
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warm-up to the bath T would give rise to a minimum in 
a{t) as observed. Likewise, a nearly adiabatic increase of 
Us would lead to the heating of the 2DES, resulting in a 
maximum in a{t) [l7| . The time-scale for equilibration 
would be determined both by heat leaks from the bath 
and by slow processes [i^ related to collective rearrange- 
ments of the 2DES to adjust to the new conditions. The 
applied | A Vg | are expected to trigger major rearrange- 
ments of the electron configuration since the correspond- 
ing shifts of the Fermi energy /^Ep ^ T 



edges the hospitality of the Aspen Center for Physics. 



13|,l2l|. This 



situation, coupled with possibly substantial changes in 
the screening of the 2DES across the MIT , is funda- 
mentally different from the cooling of the 2DES at a fixed 
Us when AT <C Ep- Charge rearrangements have been 
also invoked to interpret the nonmonotonic relaxation in 
insulating granular metals 23|. The cooling (heating) 



scenario proposed here is consistent with the data but 
other interpretations may be possible, especially since OS 
is found in many systems, ranging from manganites 24 1 
to biological systems [i^. For example, OS may be a gen- 
eral feature of systems that are far from equilibrium [26| . 

The adiabatic cooling (heating) effects should not be 
relevant for ^ Teg, i.e. for "small" excitations, since 
the final state will have a large configurational similar- 
ity with the original state due to the shortness of t^. 
Such small perturbations of a Coulor nb g lass are hence 
expected to lead to memory effects [27[, all of which 
are consistent with our results. Studies of other elec- 
tron glasses, where no OS but aging and memory effects 
were seen [5|, were done in such a perturbative regime. 
An equivalent (deep in the tw Teq limit) study in the 
2DES would be of great interest, especially across the 
MIT, but that is beyond the scope of this paper. 

Finally, we have attempted to determine the micro- 
scopic origin of the activation energy Ea ~ 57 K by 
repeating some of the measurements for different Vsub- 
In addition to changing the disorder, Vgub moves the po- 
sition of the 2D subband with respect to the bottom of 
the Si conduction band, and it also affects the splitting 
between the subbands of the 2DES [li, lli. If the ex- 



ponential process at long t results from an activation to 
Si-Si02 interface traps or to an upper subband, it can be 
estimated that the apphcd range — 5 < Vsub{y)< 

will have a significant impact on the value of Ea. We find, 
however, no change in EA^ which seems to rule out the 
above two processes as mechanisms for equilibration. 

We have demonstrated that the 2DES in Si exhibits 
aging effects, one of the hallmarks of glassy behavior, and 
identified precisely the conditions that lead to memory 
loss and nonmonotonic response. The observed complex 
dynamics of the electronic transport is strikingly similar 
to that of other systems that are far from equilibrium. 
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